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We have investigated the structural phase transitions of the transition metal oxide perovskites 
SrTiC>3, LaA103 and LaTiOa using the screened hybrid density functional of Heyd, Scuseria and 
Ernzerhof (HSE06). We show that HSE06-computed lattice parameters, octahedral tilts and ro- 
tations, as well as electronic properties, are significantly improved over semilocal functionals. We 
predict the crystal field splitting (Acf) resulting from the structural phase transition in SrTiOs 
and LaAlOs to be 3 meV and 10 meV, respectively, in excellent agreement with experimental re- 
sults. HSE06 identifies correctly LaTiC>3 in the magnetic sates as a Mott insulator. Also, it predicts 
that the GdFeC>3-type distortion in non-magnetic LaTiOs will induce a large Acf of 410 meV. 
This large crystal-field splitting associated with the large magnetic moment found in the G-type 
antiferromagnetic state suggest that LaTiOs has an induced orbital order, which is confirmed by 
the visualisation of the highest occupied orbitals. These results strongly indicate that HSE06 is 
capable of efficiently and accurately modeling perovskite oxides, and promises to efficiently capture 
the physics at their heterointerfaces. 



PACS numbers: 71.15.Mb,71.15.Ap,77.80., 77.84. 

I. INTRODUCTION 

Heterointerfaces between metal oxides like 
SrTi0 3 (STO), LaAlOa (LAO) and LaTi0 3 (LTO) 
show considerable promise as components in all-oxide 
electronics because they exhibit unique properties 
unobserved in the corresponding isolated parent com- 
pounds. In their bulk phase, SrTiOs and LaAlO-3 are 
non-magnetic wide-bandgap materials, but when assem- 
bled into superlatticcs, interesting properties such as 
high T c superconductivity, magnetism, ferroelectricity 
and colossal magnetoresistance are observed. Despite 
recent experimental and theoretical activities, the field 
of oxide interfaces remains full of surprises and unsolved 
problems. 1 ' 2 

A number of hctcrointcrface behaviors have been ex- 
plained using density functional theory (DFT), but sev- 
eral authors (see Refs. 3 and 4 and references therein) 
have shown that more elaborate methods like dynami- 
cal mean field theory 5-7 or GW-based approaches 8 are 
needed to successfully model interfaces of strongly cor- 
related electronic systems and predict the direction of 
charge transfer. While many-body techniques in combi- 
nation with DFT methods are sufficient for this task, 
they remain computationally demanding, particularly 
when it is necessary to account for the important struc- 
tural relaxations and phase transitions occurring at in- 
terfaces. 

Theory's predictive power can be enhanced by con- 
sidering the properties of the two (or more) parent ma- 
terials in the bulk phase, followed by mapping to ob- 
tain heterointerface properties. 9,10 To successfully model 
the perovskites considered herein, good band gaps and 



accurate structural/geometric properties beyond reason- 
able lattice parameters are required. For example, the 
TiOg rotation and tilt angles are likely responsible for the 
metal-insulator transition in oxide superlattices 11 and ne- 
cessitate accurate theoretical structures. 

Controlling the octahedral rotations relies on good ge- 
ometric values of TiOg and is considered by some to 
be the key to designing functional metal oxides. 12 Fur- 
ther, the TiOg rotation and/or tilts affect the electronic 
properties of both bulk metal oxides and their super- 
lattices. For example, the degeneracy of the t 2g states 
in (LAO) 5 /(LTO)„/(LAO) 5 superlattices is lifted by the 
crystal field induced by the TiOg octahedral rotation in 
LaTiOa layer. 13 Consequently, accurate determination of 
crystal field splittings (CFS) induced by the phase tran- 
sition in bulk metal oxides is also important to the Ti-3d 
orbital reconstruction at the heterointerfaces. 14 

Recently, Jalan et al. found that stress also ex- 
erts a pronounced influence on the electron mobility in 
STO, with moderate strains resulting in a greater than 
300% increased mobilities. Good theoretical estimation 
of strain in metal oxides is thus recommended and com- 
putational methods approaching the experimental strain 
are preferred. This is the case for the recent HSE06 cal- 
culations of Janotti et al., 16 which confirmed the exper- 
imental findings and showed that strain in STO affects 
seriously the electronic effective masses and the conduc- 
tivity. 

A survey of literature (see Ref. 17 and references 
therein) presented in Table I below, summarizes the main 
semilocal and hybrid functional trends for bulk SrTiOs, 
a widely-used substrate for the heterointerface super- 
lattices of interest here. Methods underestimating the 
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band gaps, like the Local Spin Density Approximation 18 
(LSD A), the generalized gradient approximation func- 
tional of Perdew, Burke and Ernzerhof 19 ' 20 (PBE), and 
the PBEsol 21 ' 22 re-parameterization for solids must be 
avoided because they lead to an overestimation of the 
two-dimensional band at the heterointerfaces, and tend 
to overestimate the octahedral rotation angle in SrTiC>3. 

In addition, LSDA and PBE also fail to describe 
the antiferromagnetic insulating character in strongly- 
correlated systems 

DFT+C/ methods can open the band gap to the ex- 
perimental value by empirical tuning of the value of the 
Coulomb on-site (U) and exchange (J) interaction pa- 
rameters, but at the cost of incorporating parameters 
that are external to theory and are material-dependent. 
While DFT+U can correctly describe the insulating char- 
acter of LaTiC>3 in the ground state, 28 geometric issues 
like the rotation angle overestimation problems in STO, 
which originate from the parent functionals LSDA and 
PBE, do not improve by adding the +U correction. De- 
spite these well-known limitations, DFT+C7 has been the 
method of choice for the oxide superlattice calculations, 
mainly because the corrected gaps approach the experi- 
mental band gaps of SrTi03, LaA103 and LaTiOs bulk 
phases with a minimal additional cost compared to reg- 
ular DFT calculations, especially for superlattices. 29,30 

Very important progress has been made in further- 
ing our understanding of metal oxides superlattices us- 
ing DFT+C/; however, one value of U will not produce 
the correct band gap for differing bulk metal oxides, 
which poses a problem for heterostructure systems like 
LAO/STO, LTO/STO, LTO/LAO, etc. 

In fact, most DFT+U calculations assume that 
the value of U appropriate for the substrate is also 
valid to describe the grown material, as well as the 
heterointerface, 11,28,31 although in some instances the U 
correction is applied to the bulk only, while the grown 
material is treated at the regular DFT level. 32 (Even 
more sophisticated methods use different values of U 
and J for the Ti-3d and La-4/ states; see for example 
Ref. 33). Another problem recently gaining notoriety is 
the fact that DFT+C/ predicts exaggerated octahedral 
distortions in the fully-relaxed superlattices, 33 which se- 
riously affects the electronic structure and the properties 
of the 2-D electron gas at the interfaces. All of these 
shortcomings indicate the need for a universally applica- 
ble functional to better describe the electronic structure 
of thin films and their interfaces. 

Electronic structure calculations on bulk SrTiC>3 using 
the global hybrid functionals B3PW, 34 B3LYP, 35,36 and 
Bl-WC 37 were performed by several groups, 24,25,38 show- 
ing improved band gaps compared to semilocal DFT. For 
example, band gap deviations from experiment were only 
6-7% with B3PW, while structural properties and order 
parameters of the antiferrodistortive (AFD) phase (see 
Table I) were also closer to experiment. Nevertheless, 
the use of global hybrids in studying heterointerfaces re- 
mains limited due to their high computational cost. 



Screened hybrid density functionals like that of Heyd, 
Scuseria, and Ernzerhof 39,40 (HSE/HSE06) and the 
middle-range screened hybrid functional of Henderson, 
Izmaylov, Scuseria and Savin (HISS) 41 are excellent can- 
didates for this task due to their accuracy and much lower 
computational cost compared to regular hybrids. 42-45 
Screened hybrid functionals yield SrTi03 structural 
properties for both phases 17,46 in very good agreement 
with experiment, 47 especially if used with large (but still 
computationally tractable) basis sets. The best results 
arise from the HSE/TZVP combination. 17 

Previous calculations on iron and transition metal 
monoxides 49 using HSE03 50 indicate that the magnetic 
moment, exchange splitting and band width in metals 
are overestimated relative to experiment. Subsequent 
tests on strongly-correlated magnetic materials have been 
carried out by Rivero et al. 23 using the short-range 
screened hybrid HSE06 and the long-range screened hy- 
brid LC-wPBE. 51 They demonstrated that both range- 
separated hybrids give a good quantitative description of 
the electronic structure of strongly-correlated systems, 
with magnetic coupling constants that are larger than 
experiment, yet improved compared to the results of the 
global hybrid, B3LYP. 

A more recent work form Philipps and Peralta in- 
volving transition metal molecular complexes suggests 
that HSE06 is a promising alternative for the evalua- 
tion of exchange couplings in extended systems, again 
because of reasonable accuracy at reduced computational 
expense. The present paper is a continuation of ear- 
lier work 17 with SrTiC>3: herein, we assess the accu- 
racy of HSE06 for calculating the crystal field splitting in 
SrTiC>3, the geometries of LaA103 in its different phases, 
as well as the structural and magnetic properties of the 
strongly-correlated system, LaTiOs. 



II. COMPUTATIONAL METHODS 

All calculations presented in this paper were per- 
formed using a development version of the Gaussian 
suite of programs, 53 with the periodic boundary condi- 
tion (PBC) 54 code used throughout. Unless otherwise 
noted, crystal structures were downloaded as CIF files 
from the ICSD. 55 (See Ref. 55 for the specific ID num- 
bers.) 

The Def2- 56 series of Gaussian basis sets were opti- 
mized for use in bulk LaTiOs and LaA103 calculations, 
following the procedure described in Ref. 17 for bulk 
SrTi0 3 . We use the notation TZVP and SZVP to differ- 
entiate these optimized PBC basis sets from the molecu- 
lar Def2-TZVP and Def2-SZVP basis sets. The function- 
als applied in this work include the LSDA 18 (SVWN5), 57 
PBE, and HSE06. LSDA and PBE calculations were uti- 
lized to assess the quality of the Gaussian basis sets via 
comparison with planewave calculations from the litera- 
ture. 

Most numerical settings in GAUSSIAN were left at the 
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TABLE I. Trends observed in the performance of DFT and post-DFT calculations applied to bulk SrTiOa regardless of the 
basis sets used (see Ref. 17). Hybrids account for both global (B3PW, 24-26 B3LYP 26 and Bl-WC. 27 ) and screened (HSE and 
HISS) hybrid functionals. 



Functional LSDA PBEsol PBE DFT+U Hybrids 



Lattice parameter Too Low Good Too Big Too Low Good 

Band gaps Too Small Too Small Too Small Good Good 

-40% -40% -40% ±6% 

Ti0 6 rotation angle Too Big Too Big Too Big Too Big Good 

Up to 8° 5° 4° 4-5° 1-2° 



default values, e.g., geometry optimization settings, in- 
tegral cut-offs, fc-point meshes and SCF convergence 
thresholds. For LaAlC>3, the reciprocal space integration 
used a 12x12x12 fc-point mesh for the cubic unit cells of 
5 atoms, while for the larger rhombohedral supercell of 
30 atoms, the default fc-point mesh of 8x8x4 was found 
to be sufficient. Because we performed optimizations of 
both atomic positions and lattice parameters, the SCF 
convergence threshold was set to tight, or 10 -8 atomic 
units. The fully relaxed structures can be obtained from 
the Cambridge Structural Database. 58 

Due to the metallic nature of LaTiC>3, a mesh of 
24x24x18 for the fc-point sampling was required in order 
to ensure stable convergence. 48 For the sake of computa- 
tional efficiency, all calculations on LTO were carried out 
using the SZVP basis set. 



III. RESULTS AND DISCUSSION 
A. SrTi0 3 

Under ambient conditions, bulk SrTiOa crystallizes in 
a cubic perovskite structure, subsequently undergoing a 
second order phase transition at T c =105 K to a tetrago- 
nal structure with slightly rotated oxygens around the 
z-axis, known as the antiferrodistortive (AFD) phase. 
The phase transition of STO is governed by two order 
parameters. The primary order parameter is the rota- 
tion angle of the Ti06 octahedra (9) that can reach up 
to 2.1° at 4.2 K. 59 The second order parameter measures 
the tetragonality of the unit cell, as defined by the ratio 
of c/a, and can be as large as 1.0009 upon cooling to 
10 K. 60 

The finer details of the conduction band splitting in- 
duced by the AFD phase transition in STO are still un- 
der debate, 61-64 with discrepancies between theory and 
experiment arising from differences in the sample dopant 
or defect concentrations: the TiOg rotation angle may 
be enhanced or reduced depending on the doping con- 
ditions and this affects directly the rotation of oxygen 
atoms in the xy-plane when STO distorts into the AFD 
phase. 

We have shown 17 in modeling an idealized STO (with 



no defects, doping, strain or surface effects) that HSE06 
provides band gaps and phase transition order parame- 
ters in excellent agreement with experiment. We found 
a c/a=1.0012 corresponding to a +0.12% tensile strain 
along the [001] direction or the z-axis. However, we did 
not address the crystal field splitting (Aqf) of the con- 
duction band minimum (CBM) induced by the phase 
transition. As is evident from the enlargement of the 
region around the calculated CBM at the Y point the 
crystal field causes the 3-fold degenerate band to split 
into a doubly-degenerate band and a single band 3 meV 
higher in energy. 

Further discussion arise from figure 1 by zooming at 
10% from Brillouin zone center along the high symmetry 
directions, [100] (A), [001] (£>), [101] (S), [110] (£) and 
[111] (A) in the body centered tetragonal Brillouin zone. 
Along A, D, and S, we observe a qualitative behavior 
very similar to that predicted in an earlier augmented 
plane wave (APW) calculation without spin-orbit (SO) 
coupling. 64 However, differences arise along £ and A di- 
rections where we observe that the lowest doubly degen- 
erate band splits into a heavy electron (he) and light 
electron (le) bands as we move farther from the Y point. 
The differences might come from the model used in the 
APW calculation by Mattheiss neglecting both the near- 
est and second-nearest neighbor oxygen-oxygen 2p inter- 
actions and the interaction between the oxygen-2s and 
the titanium t2 g orbitals in the cubic state. 

Quantitatively, the APW calculations predict 
Ap™ =20.7 meV which is much greater than the 
Aqp E06 — 3 meV we obtain in the present work. This 
is also the case for the spin-orbit splitting for cubic STO 
A$g w =83 meV, if compared with Af ( f B06 =28 meV 
computed using the VASP planewave code by Janotti et 
al.. 16 

Janotti et al. studied with HSE06 the effect of ±1% 
strain on the CBM focusing on the splittings and the elec- 
tron effective masses taking into account the spin-orbit 
effects. They also found that the rotation of TiOg oc- 
tahedra around the z-axis without including the tetrag- 
onal distortion (c > a) lead to a splitting of the doubly 
degenerate band of 4 meV with minor changes on the 
effective masses. Going back to our results in figure 1, 
one concludes that the small +0.12% tensile strain along 
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FIG. 1. (Color online) Enlargement of the region around calculated conduction band minimum (CBM) at the F point showing 
the splitting in SrTiOscaused by the AFD phase transition. Fhe high symmetry directions are those of the body centered 
tetragonal lattice [100] (£>), [001] (A), [101] (S), [110] (E), [111] (A) and Z =(00|) M =(|0|). 



the [001] direction arising from the tetragonal distortion 
might not contribute significantly to the CBM splitting 
but plays a role on the important differences in the elec- 
tron effective masses we observe along the high symme- 
try directions. For example, the curvature of the CBM 
is higher (i.e. lower effective mass or light electrons-Ze) 
along S and A directions -parallel to the strain- indicat- 
ing higher electron mobilities in those directions. This 
result is in line with the enhanced electron mobility Jan- 
otti et al. concluded for SrTiOa subject to +1% tensile 
strains along the (001) direction. However, further in- 
vestigation is still needed to draw a final conclusion since 
the electron effective masses and the order of the bands 
change upon the inclusion of SO effects. 16 

The combined HSE06 results, Ag| £06 and AgQ Eoe 
show an excellent quantitative agreement with experi- 
mental data obtained from single domain Shubnikov-de 
Haas oscillation Raman spectroscopy measurements, 63 
namely: A^=2 meV and A^q=18 meV. However, they 
do not support the picture proposed by the recent an- 
gle resolved photoemission spectroscopy (ARPES) exper- 
iment of Chang et al., 62 which indicates that the single 
band is of lower energy, with a doubly- degenerate band 25 
meV higher. This disagreement may arise from the fact 
that STO samples used in this experiment were doped 
with oxygen vacancies (Vo), which are known to intro- 
duce a shallow defect level in the band gap of STO. 65 
Finally, using HSE06 and the TZVP basis set to model 
STO relaxed in the AFD phase, we found that the crystal 
field splits the CBM at the gamma point by 3 meV, leav- 
ing it doubly degenerate in the absence of SO coupling. 



B. LaA10 3 

Bulk LaA103 undergoes a phase transition from the 
simple cubic structure to a rhombohedral-central hexag- 
onal structure with space group R3C at temperatures be- 
low 813 K. 66-68 The two order parameters for this phase 
transition are the tilt of the A106 octahedra (O), which 
can reach a maximum value of 5. 7° at 4.2K, 67 and r, a 
measure of spontaneous strain, defined as r = c/a — 



with values of up to -0.008. The calculated lattice pa- 
rameters relative to experiment for both phases of bulk 
LAO are underestimated with LSDA, overestimated with 
PBE and in-between for HSE06 (Table II). This is well- 
known behavior for LSDA and PBE (see e.g. Ref. 69 
and references therein) , and has been observed in a large 
number of semiconducting materials, as well as in our 
recent study on SrTi0 3 . 17 

In the low temperature rhombohedral phase, the LSDA 
and PBE lattice parameters tend to be higher than previ- 
ous planewave calculations. 71-73 The octahedral tilt an- 
gle for both LSDA and PBE is ca. 4°, which is 30% 
smaller than the maximum tilts measured experimen- 
tally. The SZVP basis set produces tilt angles 30% lower 
than those observed by planewave calculations 73 which 
is very similar to the behavior we observed in the case 
of STO: we attribute this to the use of localized ba- 
sis sets. 17 HSE06 yields lattice parameters for the cu- 
bic phase in excellent agreement with experiment, with 
a deviation of only 0.25%, and improved rhombohedral 
lattice parameters deviating by 0.5% compared to 1.5% 
with PBE. Improvement of the same magnitude was also 
obtained using the recently introduced variational pseudo 
self-interaction correction approach(VPSICo). 74 

The calculated spontaneous strain (r) is also signifi- 
cantly closer to experiment for HSE06 and LSDA than 
PBE. For all functionals considered here, the calculated 
lattice mismatch between cubic LAO and STO agree with 
the 2% value reported in the experiment, indicating that 
the basis set we used will not introduce any additional 
strain in LAO/STO superlattices beyond those inher- 
ently present in the experimentally measured ones. 

Turning now to the electronic properties of LAO, 
our computed band gaps are summarized in Table III. 
Comparisons with experimental band gaps are relative 
to the room temperature rhombohedral phases since, to 
our knowledge, no experimental measurements of the cu- 
bic phase band gap exist. 

In general, band gaps calculated using the LSDA and 
PBE are underestimated in average by ~30% or ~1.8 eV 
compared to experiment, with PBE doing better than 
the LSDA. This is consistent with earlier calculations 71,72 
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TABLE II. Calculated lattice parameters of LAO as well as phase transition order parameters: the octahedral tilt angle (O) and 
the spontaneous strain (r = c/a — V6). Values with a reference number come from previous experimental and/or planewave 
calculations. 





a (A) 


b(A) 


Rhombohedral 

c(A) 


60 


T 


Cubic 
a=b=c (A) 


Experiment 3 


5.365 


5.365 


13.111 


5.7 


-0.008 


3.810 b 


LSDA 


5.341, 5.290 ° 


5.339 


13.020 


4.1, 6.1 c 


-0.012 


3.760, 3.750 c 


PBE 


5.448, 5.370 d 


5.444, 5.370 cl 


13.250,13.138 d 


4.0 


-0.017 


3.830 


HSE06 


5.398 


5.393 


13.160 


3.1 


-0.011 


3.800 



a At 4.2 K from Ref 67 
b Ref 70 
c Ref 71 
d Ref 72 
c Ref 73 



TABLE III. Calculated band gaps of LaA10 3 using the SZVP 
basis set compared to experimental values. Data from this 
work are listed first, followed by data from previous planewave 
simulations and citations. 





Rhombohedral 


Cubic 




Band gap(eV) 


Direct 


Indirect 


Direct 


Experiment 


5.60 a 






LSDA 


3.75, 3.87 c 


3.25, 3.3 b 


3.46 


PBE 


3.98, 3.95 c 


3.27, 3.1 f 


3.54 


HSE06 


5.24 


4.77 


5.04 



a Ref 75 
b Ref 71 
c Ref 72 
d Ref 76 
e Ref 73 
f Ref 77 



and indicates that the basis set does not negatively im- 
pact electronic structure properties in a significant man- 
ner, a behavior we observed in the STO case as well. 17 
(Please see Table II for comparisons with previous cal- 
culations.) As expected, HSE06 provides a much better 
estimate of the band gaps, even when using the smaller 
SZVP basis sets, with deviations from experiment not 
exceeding 0.36 eV or 6% for the rhombohedral phase. 
This error is larger than the 1% deviation observed in 
the SrTiC>3 case, but is well-within the 0.3-0.4 eV range 
for non-ionic bulk semiconductors. 42 ~ 44 ' 78 

For the cubic phase, HSE06 produces an indirect band 
gap of 4.7 eV, which agrees well with the calculated gap of 
4.4 eV obtained using the screened exchange method 77 ' 79 
(sX) and planewave basis sets and the 4.61 eV gap ob- 
tained using the VPSIC method. 74 The HSE06 direct 
gap is 5.0 eV, which lies between the VPSICo' 4 value of 
4.89 eV and the GW-corrected LSDA band gap at the 
experimental lattice constants 76 of 5.68 eV. The quality 
of the cubic phase results is encouraging for future het- 
erostructure and defect state calculations because LAO 
films grown on Si substrates, using a few STO layers 
as a template, usually adopt the cubic structure. 80 The 
cubic and rhombohedral phase band structures as calcu- 



lated by HSE06 for LaA103 are shown in Figure 2. They 
look very different because the Brillouin zone sampling 
and the number of atoms in each simulation supercell 
are different, but as we will demonstrate with the pro- 
jected densities of states (PDOS), only small changes in 
the spectrum occur at the phase transition. 

Starting from the cubic phase, the band gap is indi- 
rect (R— > r) with a doubly degenerate CBM and a triply 
degenerate valence band maximum, VBM. A close look 
at our computed CBM in the rhombohedral phase also 
reveals that the phase transition induces a lift of degen- 
eracy at the CBM and a splitting at the T point with 
Aci?=10 meV. The band gap increase and Acf for LAO 
are much higher than the values we reported for SrTiOa; 
we attribute this to the fact that six of six oxygens are 
experiencing significant rigid octahedral tilts leading to 
changes in the La-O-La angles from the ideal 180° to 
173.8°, compared to a smaller rotation resulting from a 
4/6 oxygen ratio in the STO case. 

By aligning the VBM of both phases (see the right of 
Figure 2) we observe that the octahedral tilt leads to 
a shift of the CBM to higher energies, thus increasing 
the band gap by 500 meV. There is also experimental 
evidence that the crystal field caused by A106 tilts in 
LaAlOa induces splitting of the CBM; but no quantita- 
tive values were reported. 6 

By looking at the PDOS (Figure 3) we observe that 
the VBM is strictly dominated by 0-2p states, while La- 
5d states dominate the CBM. Unlike the SrTiOa case, 
no significant intermixing between the 0-2p and La-5d 
states is observed, indicating that the bonding has pre- 
dominately ionic character. This picture is not affected 
by the phase transition, and no signs of orbital inter- 
mixing is observed in the rhombohedral phase PDOS as 
well, where peaks conserve their height and shape (not 
shown) . 

In summary, by applying HSE06 to LAO we lose some 
agreement with experiment on the octahedral tilts but 
we gain better precision than semilocal functionals in the 
band gaps, lattice parameters and strain altogether. In 
other words, overall HSE06 provides a more physically 
accurate picture. 
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FIG. 2. (Color online) Band structure of LaAlOs calculated using HSE06 for the high temperature cubic phase (Left) and the 
low temperature rhombohedral phase (Right). The dashed line depicts the Fermi level at the valence band maximum (VBM). 



C. LaTi0 3 

LaTiC>3 adopts a Pbnm, GdFe03-type orthorhombic 
structure (a / i ^ c) with tilted and rotated TiOg oc- 
tahedra (see Figure 4) in the simulation supercell of 20 
atoms, with no experimental evidence of a simpler cubic 
structure. The tolerance factor of a perovskite compound 
ABO3, t, is defined as the ratio of the intrinsic sizes of 
the AO square and the BO2 square: 6 



namely, d 



.81,82 



t = 



V2(r B + r ) 



(1) 



where Ta, Tb and to are the ionic radii of each ion. In 
LaTiC>3, the tolerance factor is too small to favor the 
cubic structure due to the relative size of La +3 ion to the 
TiC>6 octahedra. 

Nevertheless, the octahedral tilts (0) stabilize the 
structure by shortening some La-0 bonds (see Figure 4- 
a) and causing Ti-Oi-Ti angles to deviate from 180°. Any 
neighboring pair the TiC>6 octahedra would tilt around 
the z-axis in opposite directions. The subsequent rota- 
tion with respect to the same axis (</>), is restricted to 
O2 oxygens forming the basal plane of the TiC>6 and oc- 
cur in the same sense. In addition, Cwik et al. showed 
that the negative sign of their orthorhombic distortion 
parameter (e = (b — a)/(b+ a)) demonstrates that the 
TiC>6 octahedra are distorted. The TiOg basal plane 
is rectangular-like instead of cubic with the O2-O2 edge 



lengths splitting into d 02 _ 02 and d °_ 02 with a ratio: 



r 2 -0 2 



Jong 

a Q 2 -Q 2 

r lshort 

a o 2 -o 2 



(2) 



r o 2 -o 2 reaches up to 1.04 at 8 K. The octahedron's 
basal plane angles also deviate from the ideal 90° (see 
Figure 4-c) leading to differences in the T1-O2 distances 



Ti-O: 



and d S Ti-b 2 anc ^ the ratio: 



TTi-Q 2 



Jong 
a Ti-Q 2 
ishort 
a Ti-0 2 



(3) 



LaTiC>3 undergoes a phase transition from a non- 
magnetic insulator to an antiferromagnetic Mott insu- 
lator at temperatures near T c =146 K without undergo- 
ing important changes in the structural parameters. To 
determine the initial geometry for LTO, we fully-relaxed 
a non-magnetic 20-atom LTO supercell with the PBE 
and HSE06 functionals, starting from several experimen- 
tal structures acquired at different temperatures. 70 ' 83 ' 84 
All structures relaxed to the same minimum. The lattice 
parameters of the fully-relaxed non-magnetic structure 
with both PBE and HSE06 are summarized in Table IV 
and compared to experimental data taken at T=155 K 
just above the phase transition temperature. 

Using HSE06, we observe excellent agreement with 
experiment as shown by the orthorhombic distortion 
parameter 85 (e) being of the right sign and magnitude. 
In addition, tq ^q 2 and rTi-o 2 ratios compare very well 
with experiment 70,83 an indication that the octahedral 
distortion is well reproduced. With PBE, e has a posi- 
tive sign which suggests that the TiOg octahedra elongate 
along b instead of a which qualitatively incorrect. The 
last behavior has been reported previously in Ref. 28 as 
a drawback of LSDA that have been overcome by adding 
the adequate U correction 86 (See the LSDA+J7 results in 
Table IV). Also, ro 2 -o 2 and Tti-o 2 ratios indicate that 
TiOg octahedra are nearly ideal and do not reproduce the 
distortion observed experimentally. The octahedral tilt 
(0) and rotation angles ($) are underestimated by 2 and 
0.45°, respectively with HSE06, but the underestimation 
is worse in the PBE case, reaching 3.3 and 2.85°, respec- 
tively. 

Turning to the electronic properties, the calculated 
projected density of states for LaTi03 in the non- 
magnetic state (Figure 5) are similar for HSE06 and PBE 
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TABLE IV. Lattice constants and orthorhombic distortion parameters as well as octahedral tilt (<d) and rotation angles (<E>) 
calculated using SZVP basis set for the fully-relaxed LaTiC>3 structures. Comparison is done with experimental data taken at 
155K from Ref. 70 near the magnetic to antiforromagnetic transition temperature Tm= 146 K. 





a (A) 


b(A) 


c(A) 


c 


ro 2 -o 2 d 


TTi-0 2 d 


©C) 


*(') 


Experiment 3, 


5.635 


5.602 


7.905 


-0.0030 


1.035 


1.011 


12.86 


9.69 


HSE06 


5.600 


5.534 


7.904 


-0.0059 


1.027 


1.013 


10.88 


9.25 


PBE 


5.595 


5.632 


7.938 


+0.0033 


1.008 


0.998 


9.54 


6.85 


LSDA b 


5.462 


5.524 


7.789 


+0.0056 










LSDA+[/ c 


5.586 


5.529 


7.89 


-0.0051 











a At 155K from Ref. 70 
b Ref 28 

c For the optimum values of (7=3.2 and J=0.9 eV identified in Ref. 28 
d Distances evaluated for the same set of reference atoms. 
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FIG. 3. (Color online) Top: Projected electronic densities 
of states (PDOS) for LaA10 3 calculated using HSE06 in the 
cubic phase. Bottom: Changes in the total density of states 
upon the phase transition form cubic (solid line) to rhombo- 
hedral (dashed line) mainly characterized by a shift of the 
conduction band states to higher energies. 



showing metallic behavior and a Fermi level lying at the 
middle of the Ti-3d band. 87 The Ti-3d band is separated 



a) b) c) 

o 




FIG. 4. (Color online) Tilted and rotated TiOg octahedra in 
LaTiOs viewed along: a) 010 direction showing the octahe- 
dral tilts (O) with respect to the z-axis b) 100 direction fea- 
turing the rectangular-like basal plane of the TiOe octahedra 
and c) 001 direction showing clearly the octahedral rotation 
($), the orthorhombic distortion and the different interatomic 
distances arising from it. Note that the lattice parameters b 
have been increased by 35% for a better illustration of basal 
plane distortion. 



from the 0-2p band by 3.95 and 2.87 eV for HSE06 and 
PBE, respectively. As with SrTiOs, the valence band 
(VB) is dominated by 0-2p states while the conduction 
band is dominated by Ti-3d states with some intermixing 
between 0-2p and Ti-3d orbitals, indicating a partially 
covalent bond. 

The band structure in the bottom of figure 5 show more 
details about the 0-2p and Ti-3d band separation where 
the t2 g triplet states 88 ' 89 can be distinguished by their 
colors. The GdFe0 3 -type distortion causes a large crystal 
field splittings (Acf) defined as the separation between 
the two first ti g triplet states 88 ' 89 at the T point and a 
subsequent smaller splitting (A' CF ) between the second 
and the third states. The amount of A.qf is strongly 
correlated to the orbital state of LTO (see the review 
paper of Mochizuki and Imada 90 and the recent paper of 
Krivenko 89 ) with small CFS inducing orbital liquid state 
(OL), 86,91 while large CFS favoring an induced orbital 
ordering (IOO). 74,92-94 

With HSE06, Acf is about 410 meV, somewhat larger 
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FIG. 5. (Color online) Top: Total and projected density of 
states for orthorhombic LaTiOa as calculated with HSE06 and 
the SZVP basis set. Only the contribution from the dominant 
states are represented. Bottom: The corresponding band 
structure of LaTiOa in the non-magnetic state. The Ti ti g 
triplet bands can be distinguished by their colors where the 
splittings caused by octahedral deformation can be clearly 
seen. 

than the experimental splittings of 240 meV 70 and 300 
meV 92 evaluated for non magnetic solutions. With PBE, 
we find a much smaller Acf of 61 meV, which is con- 
sistent with LSDA+[7 calculations ranging between 37- 
54 meV. 91 ' 95 This small CFS is a shared feature between 
semilocal functionals and might be at the origin of the un- 
settled theoretical description of the orbital state of bulk 
LaTi0 3 . 89 < 9 ° With HSE06, A' CF =190 meV which is high 
compared to the nearly degenerate doublet observed ex- 
perimentally. 70 This difference probably comes from the 
higher orthorhombic distortion we observe in our defect- 
free relaxed structures which is known experimentally to 
increase with decreasing defect concentration. 70 

Turning now to the properties of magnetic solutions; 
the experiments of Cwik et al., reveal almost no differ- 
ence in the structural properties between the magnetic 



and non-magnetic solutions close to the transition tem- 
perature. Based on that, we assume that the calculated 
structure in the magnetic state do not differ much from 
the non magnetic one near the transition temperature 
leading to A r ^ F miagn —A T ^, F gn . Thus, the total splitting 
of the t2g states in a magnetic state (A^° 5 ™) is a re- 
sult of the interplay between /\^™ ma 9 n ^ ^ e g^ruc- 
tural GdFeOa distortion and spin superexchange splitting 
(A SE ). 

The non-magnetic minimum was used for subsequent 
unrestricted spin calculations where spin flips were car- 
ried out to simulate, using the same parameters, the fer- 
romagnetic (FM) and G-type, A-type and C-type anti- 
ferromagnetic (AFM) states (see Table V and Ref. 99 for 
a description of the different magnetic orderings). Two 
reasons are beyond considering the different magnetic 
states: First, we want to test the ability of HSE06 to iden- 
tify the ground state because previous calculation with 
LDA+l/, 93 predicted the A-type AFM state to be the 
ground state, which is in disagreement with the G-type 
AFM state found experimentally. Second, those different 
magnetic orders might not be observed in bulk LTO un- 
der normal conditions, but might become relevant under 
strain or for LTO based heterostructures. Thus, details 
of electronic structure of LTO in all magnetic states will 
be examined with focus on AJ™° 9 ". 

In the ferromagnetic and antiferromagnetic spin ori- 
entations, PBE predicts LTO as metallic, whereas ex- 
periment has shown it to be a Mott insulator. 70 This is a 
well-known limitation of semilocal functionals and agrees 
very well with the conclusions from early DFT studies 
of LTO using planewave basis sets. 28 ' 87,100 ' 101 This has 
been overcome in the past by adding the Coulomb term 
correction (U) or via a GW correction. 97 ' 102 HSE06 per- 
forms much better than PBE predicting the G-type AFM 
ordering to be the ground state and a Mott insulator. 70 
LaTiOa also transforms to a Mott insulator for the other 
magnetic orders we considered, with an energy gain rela- 
tive to ground state of 21 meV for A-type AFM, 44 meV 
for C-type AFM and 795 meV for the FM states. 

The band structure of a Mott insulator is characterized 
by the optical gap, which is determined by the Hubbard 
splitting (Mgap) of the <i-band separation of the lower 
and upper Hubbard bands; this is also called the Mott 
gap. The other gap is called the charge transfer gap, 
which is the energy difference (A) between the filled p- 
bands and the unoccupied upper Hubbard 3c?-band. For 
LaTiOa, the experimental Mott gap is 0.1 eV according 
to the measurements of Arima et al. 103 while the CT gap 
is 4.5 eV. 

The nature of the band gap (direct vs. indirect), the 
shape and the width of the Mott bands as well as the 
amount of A"^ gn are discussed below: In the G-type 
AFM (Figure 6), the gap is indirect from Z — > T and 
measures 1.27 eV, separating the occupied lower Hub- 
bard band from the unoccupied higher Hubbard band, 
the Mott band is 1 eV, showing a A^ gn of 700 meV. In 
the A-type AFM ordering we observe a larger, indirect 
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TABLE V. Calculated HSE06 energy differences referenced to the G-type antiferromagnetic (AFM) ground state (AE), the 
Mott gap (Mgap), the charge transfer gap (A) and the magnetic moments (n) for the Ti atoms as obtained from the Mulliken 
population analysis of the spin density. The separation between the 0-2p states and the first Ti-3d states (A p d), the width of 
the Mott band (Em) as well as the value of the first splitting of the t2 9 at the V point (A^" 9 ") are also reported. 





Spin 


AE(meV) 


M gap (eV) 


A(eV) 


£t(a.u.) 


A pd (eV) 


E M (eV) 


A tot s (meV) 


AFM: G-Type 
Experiment 
Prev. Calcs. 






0.1-0.2 
1.6 C , 0.57 d 
0.77 c 


4.5 
5.2 C , 3.5 d 


0.57 a ,0.45 b 
0.89 c , 0.78 d 
0.68 e 


3.0 b 
3.1 c , 2.3 d 




120-300 c 
400 c , 230 d 
140 f , 54 g 


This Work 


t 1 
1 t 
















AFM G-Type 





1.27 


5.2 


0.94 


2.9 


1.0 


700 


AFM A- Type 


t t 


21 


1.39 


5.0 


0.90 


3.0 


0.6 


326 


AFM C-Type 


t 1 
t 1 


44 


0.93 


4.9 


0.92 


2.9 


1.4 


660 


FM 


t t 
t t 


795 


1.02 


•5.2 


0.92 


2.5 


1.8 


720 



a Ref 70 
b Ref 96 
c Ref 32 

d VPSIC results from ref 74 for non relaxed strucutres. 

c LDA+f/ results form ref 93 using the experimental LaTiOs structure at 8 K from Ref 70. 
f DFT+(/+GWA form ref 97. 
s LDA+C/ results form ref 98. 
h LDA+U results form ref 91. 



band gap from R — >• T of 1.39 eV, but a narrower Mott 
band (sharper and thinner) due to the reduced splitting 
between the t^g states, ca. 320 meV. In the C-type AFM 
spin orientation, the indirect gap from S — ¥ T has the 
smallest value, 0.9 eV, but A"^ s ™ is comparable to the 
G-type ordering value while the Mott band is larger by 
0.4 eV. In the FM state, we find a direct band gap of 
1 eV and AJ™" 9 ™ comparable to the G-type ordering value 
while the Mott band is larger by 0.8 eV. 



Our calculated Mott gaps, CT gaps and magnetic mo- 
ments for various magnetic orders are summarized in 
Table V and compared to experiments and previous 
calculations. The separation between the 0-2p states 
and the first Ti-3d states (A p d), the width of the Mott 
band (Em) as well as the value of the first splitting of 
the t2 g at the F point (A^" 9 ™) are also reported. In 
the G-type AFM ground state, our HSE06 value for 
A pc (=2.9 eV is in excellent agreement with the 3 eV 
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FIG. 7. (Color online) Isosurface of the highest occupied 
t2 9 state orbitals for G-type AFM LaTi03 at the P special 
k point. Orbitals show a ferro-orbital ordering alternating 
between d xz and d yz orientations between each neighboring 
TiOe octahedra. Arrows indicate the antiferromagnetic spin 
orientation on Ti sublattice. 

value reported experimentally. 96 Also, the CT gap is 
only 0.5 eV higher than the experimentally measured 
value of 4.5 eV, which constitutes an improvement over 
the previous LDA+f/ calculation. 93 However, M gap and 
li are higher than the experimentally reported values. 
This overestimation obtained with HSE06 was addressed 
in the recent study of Rivero et a/., 23 which demon- 
strated that screened hybrids do provide a good quanti- 
tative description of the electronic structure of strongly- 
correlated systems, but that the magnetic coupling con- 
stants remain overestimated 104-106 compared to exper- 
iment. This phenomenon is observed using other the- 
oretical methods, 74 ' 87,91 ' 93,97 ' 98 summarized in Table V 
and several explanations for this wide-spread divergence 
of theory and experiment have been proposed. 74 ' 90 It is 
very likely that the experimental samples are not pris- 
tine and contain defects which "metallize" them, yielding 
smaller gaps. 

Based on our HSE06 large Acf and \i values, LTO 
is expected to display an induced orbital order (IOO). 
Indeed, orbital ordering in G-type AFM LaTi03 can be 



clearly seen in Figure 7 showing the highest occupied t2 9 
state orbitals for G-type AFM LaTi03 at the T special 
k point. Orbitals show a ferro-orbital ordering oriented 
parallel to the 2-axis and oriented along the long T1-O2 
bonds. Because the long TFO2 bonds (see Figure 4-c) are 
perpendicular between each neighboring TiOg octahedra, 
the orbitals sitting along those bonds appear alternating 
between d xz and d yz and tilted due to the GdFe03 distor- 
tion. This picture confirm once again experiments 92 and 
theoretical calculations 90 suggesting that LaTi03 has an 
induced orbital order and might not be an orbital liquid. 
IV. CONCLUSIONS 

The performance of the screened hybrid HSE06 was 
tested for modeling the fully relaxed semi-conducting 
metal oxide perovskites SrTi03 and LaA103 in their 
different structural phases, and the strongly-correlated 
Mott insulator LaTi03 in the ferromagnetic and var- 
ious antiferromagnetic orderings. HSE06 gives good 
band gaps and accurate structural/geometric properties. 
It demonstrated great efficiency in handling structural 
phase transitions, and yielded accurate structural prop- 
erties and order parameters (octahedral rotations/tilts 
and strains) for all three perovskites. This is a substan- 
tial improvement over DFT+f/ - one that suggests that 
HSE06 may be as successful in treating other metal ox- 
ides and the metal-insulator transition in oxide superlat- 
tices. In addition, the crystal field splitting (Acf) of the 
ti g states resulting from the phase transitions in SrTi03, 
LaA103 and LaTi03 was evaluated for all three materi- 
als, and showed excellent agreement with experiment in 
the STO case. 63 The LTO computed A CF = 410 meV is 
rather high and indicate that LaTiOs has an induced or- 
bital order that shows clearly from our analysis of highest 
occupied X,i g state orbitals. 
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